Dynamical symmetries in Kondo tunneling through complex quantum dots 
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Kondo tunneling reveals hidden SO(n) dynamical symmetries of evenly occupied quantum dots. 
As is exemplified for an experimentally realizable triple quantum dot in parallel geometry, the 
possible values n = 3, 4, 5, 7 can be easily tuned by gate voltages. Following construction of the 
corresponding o„ algebras, scaling equations are derived and Kondo temperatures are calculated. 
The symmetry group for a magnetic field induced anisotropic Kondo tunneling is SU(2) or SO (4). 

PACS numbers: 72.10.-d, 72.15.-v, 73.63.-b 



While theoretical predictions of the Kondo effect in 
tunneling through quantum dots (QD) under strong 
Coulomb blockade conditions Q have been confirmed 
it should be born in mind that representing a real 
nanoobject by a single localized spin S=f/2 is inad- 
equate. Ubiquitous low-lying spin excitations in few- 
electron systems cannot be ignored. Even in "classical" 
planar QDs formed in GaAs/GaAlAs heterostructures 
the Kondo physics is much richer than that employed 
in analyzing the seminal experiments ||. 

The purpose of the present work is to demonstrate 
that if low-lying spin excitations are properly incorpo- 
rated, the exchange Hamiltonian of quantum dots with 
even occupation Af unveils an unusual dynamical SO(n) 
symmetry, and to suggest experiments for its elucida- 
tion. Analysis of relatively simple QD systems indicates 
the possible emergence of higher symmetries. For exam- 
ple, Kondo tunneling may be induced by external mag- 
netic field in planar QD || , since occurrence of low- lying 
triplet exciton above singlet ground state leads to an 
SO (A) symmetry. Due to Zeeman splitting, it is reduced 
to SU(2), leading to the Kondo effect in strong magnetic 
field. Similar scenario may be realized in vertical QDs 
Pj where now the Larmor (instead of the Zeeman) effect 
comes into play. Another example is a double quantum 
dot with Af = 2 which is a spin analog of a hydrogen 
molecule H2. Here, the low lying singlet/triplet manifold 
possesses the symmetry 50(4) of a "spin rotator" [^|, 

The central (and fundamental) question is then: Is 
the physics of Kondo tunneling through complex quantum 
dots intimately related with hidden SO (n) symmetries? 
The answer given below is affirmative. Moreover, these 
symmetries can be experimentally realized and the spe- 
cific value of n can be controlled by gate voltage and/or 
tunneling strength. To be concrete, the analysis is car- 
ried out below for a triple quantum dot (TQD) in a par- 
allel geometry with Af — 4 as a neutral ground state 
(see Fig. I ). It is shown to exhibit an 50(n) symmetry, 
and the relations of tunneling strengths Vi >r and gate 
voltages Vgi,V gr with the possible values n — 3,4,5,7 
and the corresponding Kondo temperatures are explicitly 
demonstrated. This example is simple enough to allow 
the construction of the corresponding o n algebras and 
solving the poor-man scaling equations for obtaining the 
Kondo temperatures. At the same time, it paves the way 



for treating more general QD structures with arbitrary 
scheme of low lying spin excitations. 

Initially, the TQD in Fig. I is treated within an 
Anderson- type model with bare level operators d„i, en- 
ergies Si, charging energies Qi and gate voltages V g i with 
i = l,c,r for left, center and right dots. The figure also 
defines inter-dot hopping (W a ) and tunneling matrix ele- 
ments (Vq) where the notation a = l,r and a — r, I is used 
ubiquitously hereafter. It is useful to shift the energies as 
ti = Si — Vgi which can be experimentally manipulated. 
Setting the Fermi energy in the leads to be Ef = 0, the 
pertinent "Kondo limit" is determined as, > e a ^> e c 
and < e a + Q a <C e c + Qc [§• The capacitive interac- 
tion between the three dots is tuned in such a way that, 
in the absence of inter-dot hopping, the neutral ground 
state has the occupation, n a = 2, n c = n a = I while 5 
electron states cost much energy and are discarded. 

Next, the isolated dot Hamiltonian is diagonalized 
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FIG. 1: Triple quantum dot in parallel geometry. 

in the Hilbert space which is a direct sum of 3 and 4 
electron states, |A) and |A), using Hubbard operators 
X ii = | 7 )( 7 | ( 7 = A, A) 0. The four particle states 
|A) = (|A/),|A,.)) exhaust the lowest part of the spec- 
trum, an octet consisting of two singlets |5/), |5 r ), and 
two triplets \Ti),\T r ). Just above it, there is a charge 
transfer exciton \ex). The corresponding energies are, 

E Sa = e c + e a + 2e a + Q- a -2W a p a , 
E Ta = e c + e a + 2e a + Qa, 

E ex = 2e l + 2e r + Q l + Q r + 2Wif3i + 2W r f3 r . (I) 

where a = W a /A a < 1 (A a = Q a + e a - e c ). Finally, 
tunneling operators in the bare Anderson Hamiltonian 
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are replaced by a product of number changing Hubbard 
operators X XA and a combination c ka = 2~ 1 / 2 (cfc (TS + 
c kad) of lead electron operators, (k — momentum, a = 
spin projection and s, d stand for source and drain). 

With these preliminaries, the starting point is a gen- 
eralized Anderson Hamiltonian describing the TQD in 
tunneling contact with the leads, 

H A = tkbc\ ah c kah + Y E 1 X^ 

kab=s,d 7— A A 

+ (££^ A 4* AA +^Y (2) 

\ AA kaa / 

with dispersion ef-b of electrons in the leads and = 
K(A|d CTa |A). The Kondo effect atT > T K is unraveled by 
employing a renormalization group (RG) procedure [0, |^| 
in which the energies E 1 are renormalized as a result of 
rescaling high-energy charge excitations. Our attention, 
though, is focused on renormalization of Eg a ,Ex a @j- 
Since the deep central level e c as well as the tunnel con- 
stants are irrelevant variables || || , the scaling equations 
are 

irdE A /d\nD = T A . (3) 

Here 2D is the conduction electron bandwidth, T\ are 
the tunneling strengths, 

r Ta = ttpo (K 2 + wi) , r 5o -a 2 r Ta , (4) 

with a a = y/1 — 2(31 an d Po being the density of states 
at ep. The scaling invariants for equations (^), 

E* A = E A (D) - 7r- x r A In(7rU/r A ), (5) 

are tuned to satisfy the initial condition E a (Dq) = E A °'. 
Equations (^) determine four scaling trajectories E\(D) 
for two singlet and two triplet states. Note that the 
level E ex is irrelevant, but admixture of the bare exci- 
ton (n a — n a = 2) in the singlet states is crucial for the 
inequality of tunneling rates IY > F$ a (cf. [|[ ||). As a 
result, the energies Er a (D) decrease with D faster than 
Es a (D), so that the trajectories £r a (-D, r<r a ) intersect 
E Sa (D,T Sa ) at certain points D^ = D { c a) . This level 
crossing may occur either before or after reaching the 
Schrieffer- Wolff (SW) limit where E A (D) ~ D and scal- 
ing terminates pi. Hidden dynamical symmetries affect 
the Kondo tunneling most effectively when the scaling 
trajectories cross near the SW boundary E A {D C ) ~ D c . 
An example of this scenario is shown in Fig. 2. As a re- 
sult, various patterns of occasional degeneracy may arise 
depending on the initial conditions ([!]), which, in turn, 
determine the pertinent SO(n) symmetry of the resulting 
spin excitations (see below). 

The above Haldane RG procedure brings us to the 
SW limit fl, where all charge degrees of freedom are 
quenched. By properly tuning the SW transformation 
e lS the effective Hamiltonian H — e %s Hjie~ %s is of the 



s — d type 0. However, unlike the conventional case H 
of doublet spin 1/2 we have here an octet A = {A;, A r } = 
{Si, Ti, S r , T r }, and the SW transformation intermixes all 
these states. To order 0(|V| 2 ), then, 

H = J2 E Ka X A ^ + J2 Jl$a ■ S + JlrP S « ' 8 
A, a a a 

+ £ J a ' S + Jlr £ B « * S + ^kbc{ ah C kah .(6) 

a a kab 

The vector operators, S a ,M Q ,B a and the permutation 
operator P manifest the dynamical symmetry of TQD. 
Their spherical components are defined via Hubbard op- 
erators connecting different states of the octet: 



9+ 

°a 


= V2(X la ° a + A° aIa ), S~-- 


- (5 a + )t, 


SI 






M+ 


= V2(X laSa — A 5aIa ), M~ 


= (M+y, 


M z a 


= -(x a « s « +X S *°"), 




B+ a 


= yftfaX 1 * 3 ' ~a a X s ^ u ), 


Ba = (B+)t 


B Z a 


= -(a- a X°^+a a X s ^). 





Here S a are spin 1 operators with projections fi a — 
la 5 a ,l o , while M a and B Q couple singlet \S a ) with 
triplet (T a fj, a \ and (T a /x a | respectively. The permutation 
P = £ Q (A S » S » + E^^) commutes with S; + S r 

and M; + M r , while S = \ Y, k y J2aa> A'a'Taa'Cka, with 

Pauli matrices f is the conduction electron spin op- 
erator. Finally, the (antiferromagnetic) coupling con- 
stants are jj = 2V^/A a , jST = a a Jj, J lr = 
ViVrY^a^a 1 ( A a = £f ~ O- After arriving at the 




FIG. 2: Scaling trajectories resulting in an SO(5) symmetry 
in the SW regime. 

s—d Hamiltonian (^) the last stage implies the solution of 
Anderson poor-man scaling equations JTcfl for extracting 
the corresponding Kondo temperatures. The discussion 
below exhausts all possible realizations of SO(n) symme- 
tries arising in TQD. 

The most symmetric case is realized when A; = A r and 
TT r = Ty,. If all four phase trajectories E\(D) intersect 
at D = D c , the symmetry of the TQD is P x 5*0(4) x 
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50(4). The operator B transforms into PM, and the 
exchange part of H (^J) reduces to 



H m t = J T ]T(1 + P)S a -s + J ST ^(1 + P)M a ■ s. 



(8) 



The vector operators M a and S a obey the commutation 
relations of 04 Lie algebra, 

\S a j , S ak \ — ^jfen^am] [M- a j , M ak ] — iejkmSam, 

[M aj ,S a k}=ie jkm M am . (9) 

(here j, k, m are Cartesian indices). Besides, S a ■ M a = 0, 
and the Casimir operator is S 2 + M 2 = 3. This justifies 
the qualification of such TQD as a double spin rotator. 
Scaling equations for Jt and Jst are, 



dji 
dlnd 



= -2 [C7i 



dji 
dlnd 



= -47'u'a, (10) 



with ji — poJ T ,j2 = PoJ ST ,d = poD. In the limit 
of complete degeneracy the system ( |io| ) is reduced to 
a single equation, dj + /d\nd — —2(j+y for j + — ji + 
]2- Its solution yields the Kondo temperature Tko = 
D exp(— l/2j + ), which is an obvious generalization of 
that derived for a QD with 50(4) symmetry and triplet 
ground state j| § §|. The net spin of the TQD is 
also 5 = 1, and the residual under-screened spin is 
5 = 1/2. If the occasional S/T symmetry is lifted, 
5 = E S (D C ) - E T {D C ) > 0, but the TQD still con- 
serves its permutation symmetry, the Kondo tempera- 
ture is not universal anymore, since the scaling of ji ter- 
minates at D ps 5 (cf. Q]). Analytic solution of Eqs. 
( ^p| ) obtains when \S\ ^> Tko, for which ji ~ ctji and 
Tr/Tko « (T KQ /S) a . The symmetry of TQD in this case 
is P x 50(3) x 50(3). 

For the symmetric configurations considered so far, the 
properties of TQD are similar to those of DQD, supple- 
mented by the permutation operation. Much richer are 
asymmetric configurations where A/ 7^ A r , Tt, ^ Dr r . 
When E Sl ~ E Tl w E Sr < E Tr (Fig. 2), the TQD 
possesses an 50(5) symmetry. The group generators of 
the 05 algebra are the " left" vectors S/ , M/ and the vec- 
tor B (with B+ = V2(X llSr - X s " 1 '), B~ = (B+)t, 
B z = ~(X°' Sr + X Sr ° 1 )), supplemented by the scalar 
operator f = i(X s - Sl - X SlS -). Thus, 

[Sij,Sik] — iejkmSim, [Mij,Mik] = iejkmSim-, 
[Bj,Sik] = iejk m B m , [Bj,Bk] — iejkmSim, 
[Mij,Sik] = ie jkm Mim, [Mij, B k ] = ifS jk , 
[Bj , f] = iMii , [f, My] = iBj , [f, 5y] = 0. (11) 

with M ; • S/ = B • Si = 0, M ; • B = 3X S ' S -, and Casimir 
operator S 2 + M 2 +B 2 + T 2 = 4. The exchange Hamil- 
tonian now reads, 



and the scaling equations are 
dji/ dlnd = -[j 2 



•if] 



dji j dlnd 



-2jih, dj 3 /dlnd = -2j 1 j 3 , (13) 



where ji = p Ju, ji = PqJu T and j 3 = p a r Ji r - From 
Eqs.([13|) the Kondo temperature is found, 



T K \ = D exp I - 



h + xljl+ji 



(14) 



Upon increasing 6 r i = Es r (D) — Et^D) the energy 
Es r is quenched, and at S r ^> Tki the symmetry 
reduces to 50(4), with Kondo temperature, Tk = 
|5r|exp{-[7i(|5 r |) +j 2 {S r )]- 1 } (cf. |). On the other 
hand, upon decreasing S r — ET r (D) — Eg r (D) the sym- 
metry P x 50(4) x 50 (4) is restored at S r < T K0 . 

Another "exotic" symmetry, namely SO (7), is realized 
when the low-lying multiplet is formed by two triplets 
Et Ut and one singlet, say Es t . In this case the 07 alge- 
bra is generated by the six vectors of the type S,M,B 
and three scalar operators describing various permuta- 
tions. Finally, an 50(3) symmetry occurs when only 
one triplet state Er a (left or right) is relevant, and the 
03 algebra is generated by S a . The dynamical symme- 
try of TQD is thereby exhausted, and summarized by 
the phase diagram in the x,y plane with x — r r /r; and 
y = A//A r depicted in Fig. 3. The central domain of di- 
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Hint = Jij'Sj • s + J^ T Mi ■ s + a r J ir B ■ s 



(12) 



FIG. 3: Phase diagram of TQD. 

mension Tko describes the fully symmetric state (||) , and 
various regimes of Kondo tunneling correspond to lines or 
segments in the {x, y} plane. Vertically hatched domain 
corresponds to TQD with singlet ground state where the 
Kondo effect is absent. Experimental test is suggested in 
Fig. 4 which illustrates the evolution of Tk, with 5 r i ~ y 
for x = 0.96 and x = 0.7 corresponding to a symmetry 
change from P x 50(4) x 50(4) to P x 50(3) x 50(3) 
and from 50(5) to 50(4), respectively. 

In similarity with planar QDs or DQDs with 50(4) 
symmetry || |], |[ H], Kondo tunneling may be induced 
by external magnetic field B also in the non-magnetic 
sector of the phase diagram of Fig. 3 close to the 50(5) 
line. In this sector 6 — Es, r — Et < 0, and the Kondo 
effect emerges when the Zeeman splitting energy E z = 
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g/is-B « 6. Due to this compensation Eti ~ Es, r and 
the spin Hamiltonian confined to this subspace has a form 
of anisotropic SU(2) Kondo Hamiltonian 



H, 



J\\R z s z + Ji (R + s- + iTs+) /2. 



(15) 



Here J,|(5) = Jf , J±{D) = ^2[(a t Jf) 2 + (a r J; r ) 2 ]. 
The vector R is defined as, 

i?+ = ^ A Q x liS -, ir = (i?+) t , (16) 

a 

i? 2 = [X 1 ' 1 ' - J2(A 2 a X 3 ^ + A a A n X s « s *)]/2, 

a 

where A t = y/2ai J\\ (£>)/ J±(D), A r = V2a r Ji r / 'Jl(-D), 
Af+A 2 . = 1, and [Rj,Rk] = ie 3 k m Rm - The operators flUf ) 
generate the algebra 03 in the spin subspace {Si, S r , Tl{\ 
specified by the Casimir operator R 2 = (3/4) LY 1 ' 1 ' + 
J2 a (AlX SaSa + A a A a X SaSii )]. The scaling equations for 
dimensionless exchange constants read, 

dj\\/dlnd= -(1l) 2 , dj ± /d\nd = -j\\jx, (17) 



yielding the Kondo temperature, 
T Kz = D exp 



1 

C 



--arctanf- 



(18) 



where C = J (2af — l)(j||) 2 + 2(a r ji r ) 2 - Another type 



the Kondo temperature 



T Kz = D exp 









-— 


arctan 


2C 


V2 





(19) 



where C = ^/(2a 2 - l){j T ) 2 . 

To conclude, the dynamical SO(n) symmetry of Kondo 
tunneling through an evenly occupied TQD is unraveled. 
It is found that the Kondo resonance with variable Tk 
arises due to strong correlations in a central well, which 
plays a role of side-coupled dot for both left and right 
wells. The hidden dynamical symmetry manifests itself, 
firstly in the very existence of the Kondo effect in QDs 
with even N ', secondly in non- universal Tk- Its depen- 
dence on the ratios x, y of the gate voltages and tunnel- 
ing rate may be observed as peculiar conductance curve 
g(x, y) at low temperature in specific Coulomb block- 
ade windows, following the curve Tk(x, y) exemplified in 
Figs. 3,4. In a singlet spin state the anisotropic Kondo 
effect can be induced in TQD by external magnetic field. 

The theory is constructed in a single-channel approxi- 
mation for lead electrons. In a split gate geometry, more 
than one tunneling channel may arise. One may antici- 
pate that the peculiar even occupation regime of complex 
QDs then transforms into conventional odd occupation 
Kondo regime. 

This research is partially supported by ISF, BSF and 
DIP funds. 
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FIG. 4: Kondo temperature 

of field induced Kondo effect is realized in symmetric 
case of 8 = Es l r — Ex t < 0. Now two components 
of a triplet, namely -Eti, r cross with two singlet states, 
and the symmetry group of the TQD is SO (4). The 
04 algebra is formed by two vectors R and PR which 
intermix the states Si, r and T^ r . The Kondo Hamiltonian 
is also anisotropic. RG procedure similar to (0) yields 
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